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Abstract. A thermal squeezed state representation of inflaton is constructed for a 
flat Friedmann- Robertson- Walker background metric and the phenomenon of particle 
creation is examined during the oscillatory phase of inflaton, in the semiclassical theory 
of gravity. An approximate solution to the semiclassical Einstein equation is obtained 
in thermal squeezed state formalism by perturbatively and is found obey the same 
power-law expansion as that of classical Einstein equation. In addition to that the 
solution shows oscillatory in nature except on a particular condition. It is also noted 
that, the coherently oscillating nonclassical inflaton, in thermal squeezed vacuum state, 
thermal squeezed state and thermal coherent state, suffer particle production and the 
created particles exhibit oscillatory behavior. The present study can account for the 
post inflation particle creation due to thermal and quantum effects of inflaton in a flat 
FRW universe. 



1. Introduction 

According to the simplest version of the inflationary scenario, the universe in the past 
expanded exponentially with time, while its energy density was dominated by the 
effective potential energy density of a scalar field, called the inflaton. Sooner or later, 
inflation terminated and the inflaton field started quasiperiodic motion with slowly 
decreasing amplitude. The universe was empty of particles after inflation and particles 
of various kinds created due to the quasiperiodic evolution of the inflaton field. The 
universe became hot again due the oscillations and decay of the created particles of 
various kinds. Form on, it can be described by the hot big bang theory. 

The standard cosmology provides reliable and tested account of the history of 
the universe from about O.Olsec after the big bang until today, some 15 billions years 
later. Despite its success, the hot big bang model left many features of the universe 
unexplained. The most important of these are horizon problem, singularity problem, 
flatness problem, homogeneity problem, structure formation problem, monopole 
problem and so on. All these problems are very difficult and defy solution within 
the standard cosmology. Most of these problems have been either completely resolved 
or considerably relaxed in the context inflationary scenario 1J. At present there are 
different versions of the inflationary scenario. The main feature of all these 



Particle creation in the oscillatory phase of inflaton 



2 



versions is known as the inflationary paradigm. Inflationary cosmology is also widely 
accepted because of its success in explaining cosmological observations 

Most of the inflationary scenarios are based on the classical gravity of the Friedmann 
equation and the scalar field equation in the Friedmann-Robertson- Walker (FRW) 
universe, assuming its validity even at the very early stage of the universe. However, 
quantum effects of matter fields and quantum fluctuations are expected to play a 
significant role in this regime, though quantum gravity effects are still negligible. 
Therefore, the proper description of a cosmological model can be studied in terms of 
the semiclassical gravity of the Friedmann equation with quantized matter fields as 
the source of gravity. The semiclassical quantum gravity seems to be a viable method 
throughout the whole non-equilibrium quantum process from the pre-inflation period 
of hot plasma in thermal equilibrium to the inflation period and finally to the matter- 
dominated period. 

Recently, the study of quantum properties of inflaton has been received much 
attention in semiclassical theory of gravity and inflationary scenarios. [Bl Ej ■ In the new 
inflation scenario jH] quantum effects of the inflaton were partially taken into account 
by using one-loop effective potential and an initial thermal condition. In the stochastic 
inflation 9\ scenario the inflaton was studied quantum mechanically by dealing with 
the phase-space quantum distribution function and the probability distribution [10J. 
The aforementioned studies show that results obtained in classical gravity are quite 
different from those in semiclassical gravity. Such studies reveal that quantum effects 
and quantum phenomena play an important role in inflation scenario and the related 
issues. Recently, it has been found that nonclassical state formalisms are quite useful 
to deal with quantum effects in cosmology [H]-[IEj, particularly squeezed sates and 
coherent state formalism of quantum optics [T5] . 

The above mentioned squeezed sates and coherent state formalisms are zero 
temperature states. There exist a thermal counterparts of coherent and squeezed sates 
and are useful to deal with finite temperature effects and quantum effects. From the 
cosmological point of view it would be more natural to consider the temperature effects 
on the background of FRW metric. Therefore this motivates the study of thermal 
squeezed states and thermal coherent states in cosmology. 

The goal of the present paper is to study quantum and finite temperature effects 
of minimally coupled massive inflaton in the FRW universe. Hence to examine the 
thermal and quantum particle creation, in the oscillatory phase, of the inflaton in 
thermal coherent and thermal squeezed state formalisms, in the semiclassical theory 
of gravity. For the present study we follow the unit system c = G = h=l. 

2. Thermal squeezed states and thermal coherent states 

The thermo field dynamic 20\ formalism can be use to get the thermal counterparts of 
coherent and squeezed states. The main feature of thermo field dynamics is the thermal 
Bogoliubov transformation that maps the theory from zero to finite temperature. One 
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can construct a thermal vacuum | 0(/3)) annihilated by thermal annihilation operators 
and can express the average value of any observable A as the expectation value in the 
thermal [20J vacuum 

Z(P)- 1 tr[pA] = (0((3)\A\0((3)), (2.1) 

where p is the distribution function, (3 = and k Boltzmann's constant, and T the 
temperature. In order to fulfill the requirement (2.1), the vacuum should belong to the 
direct space between the original Fock space by an identical copy of it denoted by a 
tilde. Therefore 

| 0{(3)) = e~ lM | 0, 0), M = -i6(p)(a)~a^ - aa), (2.2) 

where a, a 1 are the annihilation and creation operators in original Fock space and 
a, a) are the same for the tilde space, and are obeying boson commutation relations 
[a, a>] = [a, a'] — 1, the other combinations are zero. 

The density matrix approach usually gives us a convenient method for incorporating 
finite temperature effects. Hence, various definitions of thermal coherent states (tcs) can 
be summarized by giving its density matrix and it can be written for the single mode 
case as 



(he 



D\a)e- pu,a]a D{a), (2.3) 



where a is a complex number specifying the coherent state, u is the energy of the mode, 
and 

D(a) = exp [aa) — a*a). (2-4) 

The characteristic function for single mode thermal coherent state, Qt cs , is defined 
ED by 

Qtnfav*) = exp[-/(/?)| V f + ^a-na*}, (2.5) 
where n and rj* are as independent variables and, 

m = ^n- ( 2 -6) 

Similarly the density matrix for a single mode thermal squeezed states (tss) is given 
by 

p tss = D^a)S^0e^ a ' a S(0D(a), (2.7) 



where 



S(C) = exp[(^ 2 -Ca 2 )/2]^ = re w . (2.8) 



Here r is the squeezing parameter and d is the squeezing angle. 

The characteristic function of a single mode thermal squeezed state is given by 



Q tss (v,V*) = ex P[-| V | 2 ^sinh 2 rcoth^ + f{(3)j 



(2.9) 

coshrsmhr (3u , _ { 2 { , 2 x „ „ 
coth— (e ^n + e v r\ J — rja + i] a\. 
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The density matrix for a single mode thermal squeezed vacuum (tsv)is given by 

Pt sv = 5t(0e^ a+o 5(0, (2.10) 
and the characteristic function is 

QtsviViV*) =exp[— | rj | 2 ( sinh 2 r coth + /(/?)) 

V 1 J (2.11) 

cosh r sinh r /3a; / 2 4 , 2 > 
coth — I e ^ + e ^r? 



2 2 

Though the space is direct product between the original space and identical copy 
of it, the observational quantities are the expectation values of a,a\a 2 ,a^ 2 
etc. Those quantities can be computed in thermal coherent state, thermal squeezed 
state and thermal squeezed vacuum state formalisms by applying their corresponding 
characteristic function in the following relations. 
dQ 

\rj=ri*=0i 

(2.12) 

\r)=r]*=0- 



(a) 



drj* 



(a T ) = - 



drj 

Similarly the higher order expectation values of a and a) can be also evaluated using 
the same procedure of eq (2.12). 

3. Infiaton in a fiat FRW metric 

Consider a flat Friedmann-Robertson- Walker spacetime with the line element 

d s 2 = -dt 2 + R 2 (t)(dx 2 + dy 2 + dz 2 ), (3.1) 

The metric is treated as an unquantized external field. 

The minimally coupled infiaton with the gravity, for the metric (3.1), can be 
described by the Lagrangian 

L = -R 3 (ip 2 - mV) • (3-2) 

Where overdot represents a derivative with respect to time. The equation governing the 
infiaton, for the metric (3.1), can be written as 

(p + 3—0 + m 2 Lp = 0. (3.3) 
R 

One can define the momentum conjugate to cp as, it = |4. Thus, the Hamiltonian of the 
infiaton is 

" = 5ff + 5*™V. (3.4) 

Therefore, — component the energy-momentum tensor for the infiaton takes the 
following form 

^oo = y^ + mV) . (3.5) 
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Consider the minimally coupled inflaton as the source of gravity. Therefore the classical 
Einstein equation becomes 

[r) = T^' (3 - 6) 

where Too is the energy density of the inflaton, given by (3.5). In the cosmological 
context, the classical Einstein equation (3.6) means that the Hubble constant, H = j-, 
is determined by the energy density of the dynamically evolving inflaton as described 
by (3.3). 

4. Thermal and quantum particle creation 

Since there is no consistent quantum theory of gravity available, it would be meaningful 
to consider the semiclassical gravity theory to study quantum effect of matter field 
in a classical background metric. The semiclassical approach is also useful to deal 
with problems in cosmology, where quantum gravity effects are negligible. Oscillatory 
phase of inflaton is such a situation, where one can neglect the quantum gravity effects. 
Therefore the present study can be restricted in the frame work of semiclassical theory 
of gravity. In semiclassical theory the Einstein equation can be written as 

GV = 8tt(T^) , (4.1) 

where the quantum field, represented by a scalar field 0, is governed by the time- 
dependent Schrodinger equation 

^^ = . (4.2) 

. Consider quantum inflaton as the source, then the Friedmann equation, for the metric 
(3.1), in the semiclassical theory, can be written as 

(R\ 2 8?r 1 /A . , ins 

where (H^) represent the expectation value of the Hamiltonian of the inflaton in a 
quantum state under consideration. 

The inflaton can be described by the time dependent harmonic oscillator, with 
the Hamiltonian given in (3.4). To study, the semiclassical Friedmann equation, the 
expectation value the Hamiltonian (3.4) to be computed, in a quantum state under 
consideration. Therefore (3.4) becomes 

<^> = ^<- 2 > + ^<v3 2 >. (4.4) 
The eigenstates of the Hamiltonian are the Fock states 

a' (t)a(t)\n, (p,t) = n\n, (p,t) , (4.5) 

where 

a)(t) =(p(t)n - R 3 (p(t)(p . 
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As an alternative to the n representation, consider the inflaton in thermal squeezed state 
formalism. Therefor the expectation value of the Hamiltonian (4.4) in thermal squeezed 
state can be computed as follow. 

From (2.9), (2.12) and (4.6), we get 

(tt 2 ) = _ R% 2 (« 2 - e M cosh r sinh r coth ^ 

- R V 2 (a* 2 - cosh r sinh r coth (4.7) 



+ R b (p*(p 2\a\ 2 + 2 sinh 2 r coth — + 2/(/3) + 1 



and 



(y3 2 ) = — if 2 ( a 2 — e l& cosh r sinh r coth ^ 



- v?* 2 ( a* 2 - e _i * cosh r sinh r coth ^ 1 (4.S 



+ <p*(p ^2|a| 2 + 2sinh 2 rcoth^ + 2/(/5) + lJ . 

Substituting (4.7) and (4.8) in (4.4), and the apply the result in (4.3), then the 
semiclassical Friedmann equation becomes 



'i?\ 2 = 47r 
R ) ~ 3 



(<p*<p + mV» \ 2\a\ 2 + 2 sinh 2 r coth — + 2/(/3) + 1 



-(v? 2 + m 2 v? 2 ) ^a 2 - e i,? cosh r sinh r coth (4.9) 

-(y?* 2 + m 2 v?* 2 ) ^a* 2 - e^' d cosh r sinh r coth ^ 

where, y2 and </?* satisfy eq (3.3) and the Wronskian condition 

tf(t)fc*(t)<p(t)-<p*(t)<p(tj)=i. (4.10) 

The above boundary condition, fixes the normalization constants of the two independent 
solutions. 

To solve the self-consistent semiclassical Einstein equation (4.9), transform the 
solution in the following form 

<p® = -Km, (4.H) 

i?2 



thereby obtaining 



2 



-lm>r t)=0 - (4 ' 12) 

Next, focus on the oscillatory phase of the inflaton after inflation. In the parameter 
region satisfying the inequality 

o 3-R 2 3R , . 

m > w + m< ( " 3) 
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m 



—== exp(— i I a(t)dt) 



With 



a(t) 



m z — 



3 f R\ 2 3R 3 (&(t)\ 2 la(t) 



(4.14) 



(4.15) 



4 \Rj 2R 4 \a(t) J 2a(t) ' 
By applying the transform solution (4.11) in (4.9), and also using the fact a = e l5 a, we 



obtain 



R(t) 



2tt 1 

^W 2 



~ I 13^ + ^) + a 2 + m 2 



R a 



(3u 



x ^2|a| 2 + 2sinhVcoth^- + 2/(/?) + l 



x 2a 2 cos(2<5 - 2at) - cos($ - 2at) sinh(2r) coth 



(4.16) 



R 



a 



(3u 



+a I 3— + - j I 2a 2 sin(25 - 2at) + sin(tf - 2at) sinh(2r) coth - 7r 



1/3 



The next order approximation solution of the eq (4.16) can be obtained by using the 
following approximation ansatzs 



and 



Thus we get 



a (t) = m, 
R (t) = R t^ . 



(4.17) 
(4.18) 



Ri(t) 



3^mV 
1 

~2mH 2 
2 



1 



1 



2m 2 t 2 



2\a\ 2 + 2 sinh 2 r coth ^ + 2f{(3) + 1 



2a 2 cos(25 — 2mt) — cos(-# — 2mt) sinh(2r) coth 



0u 



(4.19) 



+ 



mt 2 



2a 2 sin(25 - 2mt) + sin(tf - 2mt) sinh(2r) coth ^ 



1/3 



When 26 = $ = 2mt, then (4.19) becomes 

1 



Ri(t) 



3^mV 
1 

~2m 2 t 2 



1 



2m 2 t 2 



(3u 



2\a\ 2 + 2 sinh 2 r coth '— + 2f(f3) + 1 



2a 2 - sinh(2r) coth 



f3u 



1/3 



(4.20) 
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Next, consider the particle production of the inflaton, in thermal squeezed states 
formalisms, in semiclassical theory of gravity. First, consider the Fock space which 
has a one parameter dependence on the cosmological time t. The number of particles 
at a later time t produced from the vacuum at the initial time t is given by 

N (t, t ) = (0, <p, t | N(t) | 0, <p, t ), (4.21) 

here, N(t) = o)a and its expectation value and can be calculated by using (4.6). 
Therefor, 

(N(t)} = R 6 ipip*(ip 2 } + ipip*{TT 2 } - R 3 (p(p*(frip} - R 3 ipip*(ipn}. (4.22) 
Again using (4.6) we get 

(tt 2 ) =B?(p*(p, , s 

, 4.23 

((fTt) = R 3 (p(p*. 

Therefore, substituting (4.23), in (4.22), we get 

AUMo) = &\<p(t)<p(t ) - <p(t)<p(t )\ 2 . (4.24) 

Using the approximation ansatzs (4. 17), (4. 18) and (4.24), the number of particles 
created at a later time t from the vacuum state at the initial time t in the limit mt , 
mt > 1 can be computed and is [7] given by 



N(tt)~ 1 ( R ®Y [I (z R ^ Z R ^ 



a{t h?m + (a(t)-a( to) y 



a(t) a(t ) 
(t - t ) 



(4.25) 



Am 2 4 

To compute the particle creation in thermal squeezed state, the expectation values 
of the (vr 2 ), (y? 2 ), (Try?) and (fin) in the thermal squeezed state are required. And are 
respectively obtained by using eqs (2. 9), (2. 12) and (4.6) as follow 

(n 2 } tss = -R 6 ip 2 (t ) (a 2 - e* tf cosh r sinh r coth ^ 

-R 6 ip* 2 (t ) (^a* 2 - e-' M cosh r sinh r coth 
+R 6 (p*(t )tp(t ) (2\a\ 2 + 2 sinh 2 r coth ^ + 2/(/3) + 1^ , 

{fi 2 )tss = —^(to) — ^ cosh r sinh r coth ^r^j 

-y?* 2 (t ) (a* 2 - e~ M cosh r sinh r coth ^ J 
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+<p*(t )<p(t ) (2|a| 2 + 2smh 2 rcoth^ + 2/(/3) + l) , 



(7T0) 



tss 



-R 3 (p(t )(p(t ) ( a 2 — e 119 cosh r sinh r coth ^ 



R 3 ip* (t )fi* (t ) [a* 2 - e~ M cosh r sinh r coth ^ 



+ J R 3 ( V 9(t )y?*(to) + ¥>*(*oM*o)) ( 2M 2 + 2 sinh 2 r coth ^ + 2/(/3) + 1 



(fin) tss — ~R 3 fi(to) ( p(to) \ a2 ~ e l1? cosh r sinh r coth — 



(4.26) 



R 3 ip*(t )<p*(to) I a* 2 - e~ lif cosh r sinh r coth 



+ J R 3 (0*(t o )^(^o) + fi(to)<P*(t )) ^2|a| 2 + 2 sinh 2 r coth ^ + 2/(/3) + 1^ . 

Substituting (4.26) in (4.22), we get 

1 1 1 / R(t) \ 3 



N ta8 {t,t ) 



16 a(t)a(t ) \R(to)J 



x 



R{t) o R(t ) , <r(*) <t(*o) 



-3 



+ 



x ^2|a| 2 + 2sinh 2 rcoth^ + 2/(/3) + l 

/ A(t) gfa) (T(t) a(t )\ 2 
V i2(«0) <7(f) ^ (*0) 7 



+ ^(t) 2 - a(t Q f 



(4.27) 



x 



a 2 - e i{> cosh r sinh r coth ^ | e - 2icT( *° )to 



+ I a* 2 — e ^ cosh r sinh r coth 



3 2i<x(t )t 



Which is the number of particles produced in thermal squeezed state, at a later time i 
from the initial time to- 

By using (4.17) and (4.18) the above equation (4.27) can be rewritten as follows 



tss 



i i (t-t y 

''Am 2 4 



2\a\ 2 + 2 sinh 2 r coth — + 2/(/3) + 1 



2 

-fa 2 - e i,? cosh r sinh r coth ^ ) e" 2 *" 1 ' 



(4.28) 



a* 2 - cosh r sinh r coth — | e 2imto 
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By taking a = e l5 a, eq (4.28) becomes 



tss 



^Vo(Mo) 



Pu 



2\a\ 2 + 2 sinh 2 r coth -— + 2f((3) + 1 



-2a 2 cos(25 - 2mt ) + cos($ - 2mt ) sinh(2r) coth 



Pu 



(4.29) 



where N (t,t ) is given by (4.25). 
When a = eq (4.29) leads to 

iV tes ~JV (t, t ) [2|a| 2 + 2f(p) + l-2a 2 cos(25 - 2mt )] . (4.30) 

Which is particle creation in thermal coherent state. The same result can be also 
obtained by using eqs (2.5), (2.12), (4.6), (4.17), (4.18) and (4.22). 
When r = 0, eq (4.29) becomes 



tsv 



iN (t,t ) 



2 sinh 2 r coth — + 2f(P) + 1 
2 



+ cos(-i9 — 2mt ) sinh(2r) coth 



Pu 



(4.31) 



The eq (4.31) can be also obtained by using eqs (2.11), (2.12), (4.6), (4.17), (4.18) and 
(4.22), and is the particle production due thermal squeezed vacuum state. 

When 25 = 2mt and # = 2mt , eqs (4.29), (4.30) and (4.31) respectively become 



N tss ~N (t,t ) 



2\a\ 2 + 2 sinh 2 r coth ^ + 2f(P) + 1 



Pu 



-2a 2 + sinh(2r) coth -— 



N tcs ~ N (t, t ) 2\a\ 2 + 2f(P) + 1 - 2a 2 



and 



2 sinh 2 r + 2f(P) + l + sinh(2r) coth — 

2 



N tsv ~ N (t,t ) 

When r = a = 0, then eq (4.29) take the following form 

N th ~N (t,t ) [2f(P) + l\. 
Which is the particle creation due to purely thermal effects. 



(4.32) 

(4.33) 
(4.34) 

(4.35) 



5. Conclusions 

In this paper, we studied particle production of the coherently oscillating inflaton, 
after the inflation, in thermal coherent states and thermal squeezed states formalisms, 
in the frame work of semiclassical theory of gravity. The number of particles at a later 
time t, produced from the thermal coherent state, at the initial time to, in the limit 
mto, mt > 1 calculated. It shows, the particle production depends on the coherent state 
parameter and finite temperature effects. The particle creation in thermal squeezed 
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vacuum state in the limit mt > mt > 1 is also computed, it is found that the 
particle production depending on the associated squeezing parameter and temperature. 
Similarly the number of particles produced in thermal squeezed state also computed. It 
is observed that, when r = 0, the result agree with the number of particles produced in 
the thermal coherent state and when a — 0, the result equal to the number of particles 
created in thermal squeezed vacuum state. 

The approximate leading solution obtained for the Einstein equation, in the thermal 
squeezed sates shows oscillatory behavior except when the condition, 25 = $ = 2mt, 
satisfies. Though both classical and quantum inflaton in the oscillatory phase of the 
inflaton lead the same power law expansion, the correction to the expansion does not 
show any oscillatory behavior in semiclassical gravity in contrast to the oscillatory 
behavior seen in classical gravity only when 25 = d = 2mt. It is also noted that, 
the coherently oscillating inflaton, in thermal squeezed vacuum, thermal squeezed and 
thermal coherent states representation, suffer particle creation and created particle 
exhibit oscillations. The oscillation of the created particles is necessary to preheat 
the universe to hot again after the inflation. The present study can account for the post 
inflation particle creation due to thermal and quantum effects of inflaton in a flat FRW 
universe. Since the created particle oscillate, we hope that this kind of study can light 
on preheating issues of post inflationary scenario. 
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